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Abstract 

The semi-empirical macroscopic-microscopic mass formula is further improved by considering 
some residual corrections. The rms deviation from 2149 known nuclear masses is significantly 
reduced to 336 keV, even lower than that achieved with the best of the Duflo-Zuker models. The 
a-decay energies of super-heavy nuclei, the Garvey-Kelson relations and the isobaric multiplet mass 
equation (IMME) can be reproduced remarkably well with the model, and the predictive power of 
the mass model is good. With a systematic study of 17 global nuclear mass models, we find that 
the quadratic form of the IMME is closely related to the accuracy of nuclear mass calculations 
when the Garvey-Kelson relations are reproduced reasonably well. Fulfilling both the IMME and 
the Garvey-Kelson relations seems to be two necessary conditions to improve the quality of the 
model prediction. Furthermore, the a-decay energies of super-heavy nuclei should be used as an 
additional constraint on the global nuclear mass models. 
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I. INTRODUCTION 



The precise calculation for nuclear masses is of great importance for not only nuclear 
physics but also nuclear astrophysics. In nuclear physics, the synthesis of super-heavy nu- 
clei [1-3] and nuclear symmetry energy [4, 5] attracted much attention in recent years. 
Accurate predictions for the shell corrections and a-decay energies of super-heavy nuclei 
are urgently required for the synthesis of new super-heavy nuclei. Simultaneously, nuclear 
masses include important information on nuclear symmetry energy [6, 7] which closely re- 
lates to the structure of nuclei and the dynamical process of nuclear reactions. In addition, 
the r-process abundances in nuclear astrophysics depend strongly on the nuclear masses of 
extremely neutron-rich nuclei and the shell-structure data of neutron magic nuclei. Some 
nuclear mass models have been developed with rms deviations of about several hundreds 
to a thousand keV with respect to all known nuclear masses. The most successful global 
mass models include: 1) various macroscopic-microscopic mass models such as the finite 
range droplet model (FRDM) [8], the extended Bethe-Weizsacker (BW2) formula [9] and 
the Weizsacker-Skyrme (WS) mass models proposed very recently by Wang et al. [10, 11]; 2) 
various microscopic mass models based on the mean-field concept such as the non-relativistic 
Hartree-Fock-Bogoliubov (HFB) approach with the Skyrme energy-density functional [12] 
or the Gogney forces and the relativistic mean-field (RMF) model [13]; and 3) the Duflo- 
Zuker (DZ) mass models [14, 15] which are successful for describing known masses, with 
an accuracy at the level of 360 keV. For known masses, the predictions from these global 
mass models are close to each other. However, for neutron drip line nuclei and super-heavy 
nuclei, the differences of the calculated masses from these models are quite large. It is 
therefore necessary to perform a systematic study of these models to check the reliability 
for extrapolations. 

Besides the global mass models, some local relations and equations such as the isobaric 
multiplet mass equation (IMME) [16, 17], the Garvey-Kelson (GK) relations [18, 19] and 
the residual proton-neutron interactions [20, 21] are used to analyze the isospin-symmetry 
breaking effects and the consistency of nuclear mass predictions. It is thought that the 
binding energies of extremely proton-rich nuclei can be predicted at the level of 100 ~ 200 
keV based on the IMME [16]. As the IMME is a basic prediction leading from the isospin 
concept, testing the validity of this equation is of fundamental importance. On the other 
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hand, it is shown that the GK relations which are obtained from an independent-particle 
picture and the charge-symmetry of nuclear force, are satisfied to a very high level of accuracy 
by known masses, with an rms deviation of about 100 keV [19]. It is of crucial importance 
to incorporate these local relations and equations in the global nuclear mass models for 
exploring the missing physics and the constraints in the models. 

In our previous works [10, 11], we proposed a semi-empirical nuclear mass formula based 
on the macroscopic-microscopic method together with the Skyrme energy- density functional, 
with an rms deviation of 441 keV for 2149 known masses. The isospin and mass dependence 
of the model parameters such as the symmetry energy coefficient and the constraint be- 
tween mirror nuclei from the isospin symmetry in nuclear physics play an important role for 
improving the accuracy of the mass calculation. In this work, the formula will be further 
improved by considering some residual corrections of nuclei. With some empirical consider- 
ations for these residual corrections, we find both the precision and the consistency of the 
mass formula can be significantly improved. The paper is organized as follows: In Sec. II, 
the mass formula and some corrections are briefiy introduced. In Sec. Ill, calculated nuclear 
masses from the proposed model are presented and compared to the results from other mass 
models. In Sec. IV, a number of global nuclear mass models are systematically tested by 
using the Garvey-Kelson relations and the IMME. Finally, a summary is given in Sec. V. 

II. THE NUCLEAR MASS MODEL AND SOME CORRECTIONS 

In this mass formula, the total energy of a nucleus is written as a sum of the liquid-drop 
energy, the Strutinsky shell correction and the residual correction Ares, 

E{A, Z, 13) = Eld(A, ^) n (1 + ^^f^l) + ^' + (1) 

k>2 

The liquid-drop energy of a spherical nucleus E'ld Z) is described by a modified Bethe- 
Weizsacker mass formula, 

Eld(^, Z) = a^A + asA^I^ + Ec + a.y^l'^A + a^^^^A^^/^dnp + Aw (2) 

with the isospin asymmetry I = {N — Z)/A, the Coulomb energy, 

Ec = (1 - 0.76Z-'/') (3) 
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FIG. 1: (Color online) (a) Positions of some known doubly-magic nuclei (squares). The stars denote 
the predicted tetrahedral doubly-magic nuclei in Refs.[23, 24]. (b) Shell corrections of nuclei with 
WS [10]. The solid line shows the position of iV = 1.37Z + 13.5. 



and the symmetry energy coefficient, 
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The Cpair term empirically describes the pairing effect (see Ref. [10] for details). The terms 
with hk describe the contribution of nuclear deformation (including /32, Pa and Pq) to the 
macroscopic energy. Mass dependence of the curvatures hk in Eq.(l) is written as [10], 

hk={^^9.Ay'+{^^\.A~^/\ (5) 

according to the Skyrme energy-density functional, which greatly reduces the computation 
time for the calculation of deformed nuclei. 

In addition, we propose an empirical correction term Aw in i^LD to consider the Wigner- 
like effect of heavy nuclei [22] which is due to the approximate symmetry between valance- 
proton and valance-neutron. In [10, 11], we found that some heavy doubly-magic nuclei 
surprisingly lie along a straight line = 1.37Z + 13.5 (see Fig.l), which is called shell 
stability line of heavy nuclei. Furthermore, some predicted tetrahedral doubly-magic nuclei 
[23, 24] [stars in Fig. 1(a)] also lie near the line. The shell corrections of nuclei obtained 
with WS [10] are shown in Fig. 1(b), and they are approximately symmetric along the shell 
stability line, i.e. the shell correction of doubly-magic nucleus ^^^Pb is close to that of 
^^^Sn. In other words, the change of the shell energy by adding several neutrons to the 



doubly-magic nucleus ^°^Pb is close to that by adding several protons to ^osp^ because 
of the charge-independence of nuclear force and the fact that the Coulomb interaction is 
weak relative to the strong interaction. To consider the symmetry and interaction between 
valance-proton and valance-neutron, we phenomenologically write the Wigner-like term as, 

AH. = c^(el^l-e-l^'l). (6) 

Where, rj = ^"^^^ denotes the effective distance of a nucleus to the shell stability line. The 
values of AZ, AN and are uniquely determined by a rectangle [see the green rectangle 
in Fig. 1(a)] based on the position of a nucleus under consideration and the shell stability 
line. Similar to the known Wigner effect causing a cusp for light nuclei near the N = Z line, 
the Wigner-like term A^ causes a cusp for heavy nuclei limited to the close vicinity of the 
shell stability line, which reduces the rms deviation with respect to masses by about 5%. 
The microscopic shell correction 

AE = c,E,^ + 1/1^:^ (7) 

is obtained with the traditional Strutinsky procedure [25] by setting the order p = 6 of the 
Gauss-Hermite polynomials and the smoothing parameter 7 = 1.2/ia;o with hujQ = 41^4"^/^ 
MeV. Esh and Ei.^ denote the shell energy of a nucleus and of its mirror nucleus, respec- 
tively. The |/| term in AE is to take into account the mirror nuclei constraint [11] from the 
isospin symmetry, with which the accuracy of the mass model can be significantly improved 
(by 10%). The single-particle levels are obtained under an axially deformed Woods-Saxon 
potential [26]. Simultaneously, the isospin-dependent spin-orbit strength is adopted based 
on the Skyrme energy-density functional, 

A = Ao(l + f) (8) 

with Ni = Z for protons and Ni = N for neutrons, which strongly affects the shell structure 
of neutron-rich nuclei and super-heavy nuclei. 

Some other residual corrections caused by the microscopic shell effect are empirically 
written as a sum of three terms. 

Ares = Am + Ap + At. (9) 

The first term Am is to further consider mirror nuclei effect. In our previous work [11], 
we found that the shell correction difference \AE — AE'\ between a nucleus and its mirror 
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nucleus is small due to the isospin symmetry in nuclear physics. Here, we empirically write 

AM = C2{S + S^)il-PA), (10) 

with 

Here, T = [ [ denotes the isospin of a nucleus, and we set S" = for nuclei with N = Z. 
Si is the corresponding value of a neighboring isobaric nucleus (A, Zi) with Zi = Z + 1 for 
N < Z cases and Z — 1 for N > Z cases. The Am term effectively describes the residual 
mirror nuclei and isospin effect in nuclei and significantly improves the accuracy of mass 
calculation by about 50 keV. 

The second term Ap further considers the residual pairing corrections of nuclei which may 
be phenomenologically given by the second differences of the masses. That is, in addition 
to the apair term in E'ld for describing the pairing correction, we further consider the term 
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2 dA^ 

with a step size of AA = 2 in the calculations. This term improves the smoothness of the 
mass surface assumed in the Garvey-Kelson relations [18], and simultaneously reduces the 
rms deviation by about 4% and 6% with respect to the masses and the neutron separation 
energies, respectively. In addition, for describing the influence of triaxial (or tetrahedral) 
deformation [23, 24] effectively, we proposed a parameterized formula in our previous work 
[11]. Here the formula is slightly extended to include the contribution of the protons, = 
-0.7 [cos (27rg) + cos (27rf ) + cos {2tt§)] A~^/^ MeV. 

III. CALCULATED MASSES FROM THE MODEL 

Based on the 2149 (A^ and Z > 8) measured nuclear masses AME2003 [27], we obtain 
the optimal model parameters which is labelled as WS3 and listed in Table I. We show in 
Table II the rms deviations a{M) between the 2149 experimental masses and predictions 
of some models (in keV). cr(S'„) denotes the rms deviation to the 1988 measured neutron 
separation energies Sn- ct^Qq) in Table II denotes the rms deviation with respect to the 
a-decay energies of 46 super-heavy nuclei {Z > 106) [11]. cr(GK) presents the rms deviation 
with respect to the GK mass estimates and cr(6) denotes the rms deviation from the formula 
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TABLE I: Model parameters of the mass formula WS3. 



Quantity 


Value 


Quantity 


Value 


(MeV) 


-15.5485 


91 


0.01037 


as (MeV) 


17.4663 


92 


-0.5071 


a, (MeV) 


0.7128 


Vo (MeV) 


-45.2092 


Csym(MeV) 


29.1174 


ro (fm) 


1.3848 




1.3437 


a (fm) 


0.7617 


e 


1.1865 


Ao 


26.6744 


apair(MeV) 


-6.2299 


Cl 


0.6454 


Cw (MeV) 


1.0490 


C2 (MeV-^) 


1.6179 



TABLE II: rms a deviations between data AME2003 [27] and predictions of some models (in keV). 
The line a{M) refers to all the 2149 measured masses, the line <t(S'„) to the 1988 measured neutron 
separation energies Sn, the line (j{Qa) to the a-decay energies of 46 super-heavy nuclei. a{GK) 
presents the rms deviation with respect to the GK mass estimates. a{h) denotes the rms deviation 
from the formula [16] h^^ = 0.710A^/'^ — 0.946 for the h coefficients of 179 neutron drip line nuclei. 



ETF2[28] BW2[9] FRDM[8] HFB-17[12] DZ10[14] WS[10] WS*[11] DZ31[14] DZ28[15] WS3 



a{M) 


3789 


1594 


656 


581 


561 


516 


441 


362 


360 


336 


cr{Sn) 


1300 


586 


399 


506 


342 


346 


332 


299 


306 


286 


cr{Qa) 


557 


1233 


566 




916 


284 


263 


1052 


936 


248 


a{GK) 


24 


129 


337 


496 


100 


151 


165 


115 


133 


131 


a{b) 


2734 


759 






870 


283 


450 


630 




449 



[16] 6fit = 0.710A^/^ — 0.946 for the b coefficients of 179 neutron drip line nuclei, which will 
be discussed in the next section. Here, ETF2 denotes the extended Thomas- Fermi approach 
including all terms up to second order in the spatial derivatives [28] together with the Skyrme 
force SkP [29], in which the deformations and shell corrections of nuclei are not involved. 
BW2 means the extended Bethe-Weizsacker formula [9] in which the shell corrections of 
nuclei are described as a function of valence-nucleon number assuming the magic numbers 
being 2, 8, 20, 28, 50, 82, 126, 184. FRDM and HFB-17 denote the finite-range droplet model 
[8] and the latest Hartree-Fock-Bogoliubov (HFB) model with the improved Skyrme energy- 
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FIG. 2: (Color online) (a) Deviations of calculated nuclear masses from the experimental data as 
a function of neutron number. The squares and the solid circles denote the results of FRDM and 
WS3, respectively, (b) The same as (a), but with the results of DZ28 model for comparison. 

density functional [12], respectively. DZIO, DZ28 and DZ31 denote the Duflo-Zuker mass 
models with 10, 28 and 31 parameters, respectively [14, 15]. 

The rms deviation from the 2149 masses with WS3 is remarkably reduced to 336 keV, 
much smaller than the results from the finite range droplet model (FRDM) and the latest 
HFB-17 calculations, even lower than that achieved with the best of the Duflo-Zuker models. 
In Fig. 2, we show the deviations from the experimental masses as a function of neutron 
number. The squares and the solid circles in Fig. 2 (a) denote the results from the FRDM and 
WS3, respectively. The accuracy of mass calculation especially for light nuclei is significantly 
improved in the WS3 model. The open circles in Fig. 2(b) denote the results from the 
DZ28 model. Compared to the DZ28 model, the rms deviations from masses and neutron 
separation energies are reduced by 7%. To further test the extrapolation of our mass model, 
we also calculate the average deviation from four very recently measured masses for ^^Ge, 
^^As, ^''Se, and ^^Kr [30]. Our result is 118 keV, much smaller than the results from DZ28 
(321 keV) and FRDM (680 keV). 

It is known that the a-decay energies of super-heavy nuclei [1] have been measured with 
a high precision, which provides us with useful data for testing our mass model. The rms 
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FIG. 3: (Color online) a-decay energies of some super-heavy nuclei calculated with WS3 (solid 
curve) and DZ28 (open circles). The squares denote the experimental data taken from [11] (and 
references therein). 

deviation (T[Qa) with respect to the a-decay energies of 46 super-heavy nuclei from the 
WS3 model is remarkably reduced to 248 keV. The difference between the calculated a- 
decay energies of the super- heavy nuclei and the measured data is within AQa = ±0.4 MeV 
for odd-Z nuclei. The corresponding results from Sobiczewski [31] and FRDM are ±0.8 and 
±1.0 MeV, respectively. Here, we would like to emphasize that the experimental data of 
a-decay energies are not involved in the fit for the model parameters. We also note that the 
rms deviations (T{Qa) with the extended Bethe-Weizsacker (BW2) formula and the DZ mass 
models are large, which might be due to the difficulty in determining the right shell closures 
after ^°^Pb in these models where the shell corrections of nuclei are described as a function of 
valence-nucleon number. Fig. 3 shows the calculated a-decay energies of some super-heavy 
nuclei with the WS3 (solid curve) and DZ28 models (open circles), respectively. The results 
from the WS3 model look much better than those from the DZ28 model, although both 
models give similar rms deviation for known masses. 

In Fig. 4 we show the deviations of the calculated total energies of nuclei at their ground 
state with WS3 from the results of (a) FRDM, (b) DZ28, (c) WS* and (d) HFB-17, respec- 
tively. The shades denote the region with deviations smaller than 2 MeV. One sees that the 
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FIG. 4: (Color online) Deviations of the calculated nuclear total energies at ground state in this 
work from the results of (a) FRDM, (b) DZ28, (c) WS* and (d) HFB-17, respectively. The 
calculated masses with FRDM, DZ28, WS* and HFB-17 are taken from [8], [15], [11] and [12], 
respectively. The shades denote the region with deviations smaller than 2 MeV. 

results from WS3, FRDM, DZ28 and WS* are relatively close to each other, and the results 
from HFB-17 model deviate obviously from WS3 calculations for extremely neutron-rich 
nuclei. We also note that the binding energies of super-heavy nuclei around ~ 170 and 
Z ~ 100 — 110 with DZ28 are significantly larger by about 4 MeV than those with WS3, 
which leads to the deviations in the calculation of the a-decay energies (see Fig. 3) between 
the two models. Because the predictions from different models towards the neutron-drip 
line tend to diverge [32], it is necessary to analyze the systematic errors from the various 
mass models for nuclei near the neutron drip line. In Fig. 5 we show the rms deviations 
with respect to the masses as a function of average neutron-separation energy {Sri) of nuclei 
based on the masses from different mass models. The shades present the average standard 
deviation errors of the measured masses [27]. One sees that the rms deviations from the 
FRDM, HFB-17, DZ28 and WS3 increase systematically for nuclei approaching the neutron 
drip line (i.e., the neutron separation energy Sn ~ 0), which indicates that the systematic 
errors for the neutron drip line nuclei from these mass models expect to be much larger 
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FIG. 5: (Color online) rms deviation with respect to the masses as a function of average neutron- 
separation energy of nuclei. The shades present the average standard deviation errors of the 
measured masses [27]. 

than those for stable nuclei. The rms deviation cr(M) from the FRDM and HFB-17 are 
about 1.4 and 1.0 MeV for nuclei with Sn ~ 1 MeV, respectively. The corresponding results 
from WS3 and DZ28 are about 0.8 and 0.9 MeV, respectively. The experimental errors are 
about 0.4 MeV for the masses of nuclei near the neutron drip line, and the tendency of the 
errors for both experimental data and the model calculations is similar with decreasing of 
Sn of nuclei. For nuclei approaching the proton drip line (with large value for Sn), the rms 
deviations from DZ28 and WS3 are bout 0.45 MeV, and the result from the FRDM reaches 
about 1 MeV. Fig. 5 indicates that the model WS3 is relatively reliable for describing the 
masses of nuclei far from the /3-stability line. 

Furthermore, to check the predictive power of various models, we follow Ref.[32] and 
study the quantity Y = '^^'^^^^^2003)^^^^ ' ^ measure of the predictive power of the mass 
formulas fitted to the 1995 data. Here, (t(1995) and cr(2003) denote the rms deviation to 
the masses of the AME1995 [33] and AME2003 [27], respectively. Fig. 6 shows the values 
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FIG. 6: (Color online) A comparison of the predictive power of the various models. The vertical 
axis denotes the quantity Y = °"(^003)^^o-^(^i995) ^ ^ measure of the predictive power of the mass 
formulas fitted to the 1995 data. Here, (j(1995) and (t(2003) denote the rms deviation to the masses 
of the AME1995 [33] and AME2003 [27], respectively. See Table II for more information. 




FIG. 7: (Color online) A comparison of rms deviation with respect to the masses from AME2003 
with the various models. The bars with shades are for the 2149 measured masses and the gray bars 
denote the rms deviations to the evaluated masses (marked by 7^) in AME2003. For the IIFB-17 
model, only the masses of nuclei with Z < 110 are involved in the calculation. 
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of Y with various mass model listed in Table II. Large positive value for Y indicates the 
deterioration for prediction of "new" masses in 2003 compilation. One can see from Fig. 6 
that the predictive power of FRDM, HFB-17 and WS3 seems to be relatively better than 
those of DZIO [14] and ETF2. In Fig. 7, we show a comparison of the rms deviation with 
respect to the masses from AME2003 with the various models. The bars with shades denote 
the rms deviation to the 2149 measured masses and the gray bars denote that to about 900 
evaluated masses (marked by in AME2003. One sees that the predicted masses from WS3 
are relatively close to the evaluated masses in AME2003, with an rms deviation of 0.8 MeV 
(the evaluated masses in AME2003 are not involved in the fit for the model parameters). 
The corresponding rms deviations from the three models BW2, DZIO and DZ28 in which 
one needs to determine the magic numbers in the calculations are obviously large than the 
others, like the large rms deviations with respect to the a-decay energies of super-heavy 
nuclei from these three models [see the values of (j{Qa) in Table II]. Figs. 2, 3, 5 and 6 
demonstrate that the proposed model in this work has a relatively high accuracy and good 
predictive power for a description of the masses of the drip line nuclei and the super-heavy 
nuclei. 

IV. CONSTRAINTS ON MASS MODELS FROM LOCAL RELATIONS 

In this section, we perform a systematic test of the global nuclear mass models by using 
the Garvey-Kelson relations and the isobaric multiplet mass equation. We mainly compare 
the results from the mass models mentioned previously. We first check the Garvey-Kelson 
(GK) relations based on the obtained mass tables from these mass models. It is known that 
the GK relation can not be used to cross the N = Z line. Therefore, we suitably select five 
of the neighbors in GK relation for each nucleus to avoid crossing the N = Z line. We use 
the following GK relations: 



Mgk{N, Z) 



M{N + 1,Z) + M{N, Z -1) + M{N 



+ 



2,Z - 2) 



M{N + 1, Z - 2) - M{N + 2,Z -1) 



ioT N > Z cases 



(13) 



and 



Mgk{N, Z) 



M{N -1,Z) + M{N, Z + 1) + M{N 



2,Z + 2) 



M{N - 1, Z + 2) - M{N -2,Z + 1) 



ioT N < Z cases 



(14) 
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Here, M{N, Z) denotes the mass excess of a nucleus with neutron number and charge 
number Z. cr(GK) in Table II presents the rms deviation of the calculated masses by 
the indicated mass models from the GK mass estimates according to Eqs.(13) and (14). 
It is known that the ETF approach or the Hartree-Fock-Bogoliubov (HFB) approach [34] 
based on the traditional Skyrme force without further improvements describe the nuclear 
masses rather poorly, although the GK relations are reproduced very well. With a series 
of improvements on the force, the latest HFB-17 model [12] can reproduce the nuclear 
masses with an accuracy of 581 keV. However the deviation of the calculated masses from 
the GK mass estimates becomes relatively large as mentioned in [18]. This indicates that 
more constraints are still required for further testing the reliability of the mass models, in 
addition to the GK relations. In other words, fulfilling the GK relations is a necessary but 
not sufficient condition to improve the quality of the model prediction. 

At the same time, we study the isobaric multiplet mass equation (IMME) which is widely 
used to study the isospin-symmetry breaking effect in nuclei. In the IMME, the binding 
energy {BE) of a nucleus is expressed as 

BE{T,T,) = a + bT, + cTl (15) 

with = (A^ — Z)/2. It is known that the IMME is valid if the Coulomb interaction and 
charge-nonsymmetric parts of the nucleon-nucleon interaction are weak relative to the strong 
interaction, and it is generally assumed that the quadratic term in the equation is adequate 
[4, 16, 17]. The b coefficients for light nuclei can be extracted from the measured binding 
energies of pairs of mirror nuclei [16], 

^ BEjT = - BE{T = -T,) 

2T ^ ' 

with the isospin T = |A^ — Z|/2. Based on the liquid-drop model, one can roughly estimate 
the value of the coefficient b ~ a^A^^^. Here, denotes the Coulomb energy coefficient in 
the liquid-drop model. In [16], Ormand proposed a formula 

6fit = 0.710^2/3 _ 0.946 (17) 

for the b coefficients by fitting the experimental data oi A < 60. In Fig. 8 (a), we show the 
extracted b coefficients of nuclei in the IMME from mirror nuclei with A <To according to 
Eq.(16). Here, the experimental binding energies of mirror nuclei are taken from AME2003 
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FIG. 8: (Color online) (a) b coefficients of nuclei in the IMME as a function of mass number. 
The solid squares and the solid curve denote the extracted results from mirror nuclei with A <75 
according to Eq.(16) and those from the empirical formula [16] 6fit = 0.710^4^/^ — 0.946, respectively. 
The solid (open) circles denote the extracted results by fitting the measured binding energies with 
parabolas for a series of isobaric even-j4 (odd- A) nuclei. See text for details, (b) The same as (a), 
but with the results from some mass models for comparison. The crosses, the open circles and 
the filled circles denote the results for 179 neutron drip line nuclei from ETF2, DZIO and WS3 
calculations, respectively. 

[27] and the four very recently measured masses for ^'^Ge, ^^As, ^^Se, and ^^Kr in Ref. 
[30]. It seems to be that the b coefficients are roughly constant for a given A. The solid 
curve denotes the results of the 6fit from Eq.(17). For intermediate and heavy nuclei, the 
b coefficients could be extracted by fitting the measured binding energies of a series of 
isobaric nuclei with parabolas. Here, we remove the shell corrections [11] and the Wigner 
energies (about 47]/] MeV) [35] of nuclei from the measured binding energies of nuclei before 
performing a parabolic fitting in order to eliminate the oscillations in the results. The solid 
(open) circles denote the extracted results for a series of isobaric even-A (odd-A) nuclei, 
with 95% confidence interval given by the error bars. One can see that the trend of the 
b coefficients with A can be reasonably well described by the empirical formula bf^, which 
indicates that the IMME is valid for both light and heavy nuclei in general. The residual 
oscillations in the b coefficients might come from the deformations of nuclei. 
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According to the macroscopic-microscopic mass model, the binding energy of a nucleus 
can be expressed as 

BE{A, T,) =a + bT, + dT^ + + B^,^ + Bdcf + ... (18) 

Here, the mass-dependent coefficients a, b and c can be directly obtained from the liquid- 
drop formula, and b ~ OcA^/^. B^h, B„ig and i?dcf denote the shell correction, the Wigner 
energy and the deformation energy of a nucleus, respectively. From Eqs.(15) ~ (18) and 
Fig. 8(a), we learn that: 1) the shell correction difference between a nucleus and its mirror 
nucleus is small due to the isospin symmetry in nuclear physics as mentioned in our previous 
work [11]; 2) the Wigner energy of a nucleus equals generally to that of its mirror nucleus 
since it is usually expressed as a function of |/|; and 3) the deformation energy of a nucleus 
might be close to the value of its mirror nucleus. The mirror nuclei effects from the isospin 
symmetry could give a strong constraint on nuclear mass models. 

The crosses, the open circles and the filled circles in Fig. 8(b) denote the calculated 
results with Eq.(16) for 179 neutron drip line nuclei {Z > 8 and N < 200) from ETF2, 
DZIO and WS3 calculations, respectively. The results from WS3 are relatively close to 
the values of fegt (with an rms deviation of 449 keV), which indicates that the IMME with 
quadratic form is generally satisfied in this model. For extremely neutron-rich nuclei, the 
results (crosses) from the traditional Skyrme EDF obviously deviate from the values of bta- 
From Fig. 9(a), we note that the change of fegt — b with the isospin T for the 179 neutron 
drip line nuclei is quite different with different models. For the Skyrme EDF calculations 
the difference b^t — b increases soon with the isospin T (we draw the same conclusion with 
Skyrme HFB codes HFBRAD [36] and HFBTHO [37]). Without special consideration for 
the isospin-symmetry restoration in the traditional Skyrme EDF, the calculated surface 
diffuseness of protons becomes abnormally large in extremely proton-rich nuclei to reduce 
the huge Coulomb repulsion, which causes the large deviations from feat. These investigations 
show that the exploration of the IMME in the neutron-rich nuclei might provide us with 
very useful information for testing the mass models. 

Closely relating to the IMME, the Coulomb energies (direct term) of nuclei are in- 
vestigated simultaneously. Fig. 9(b) shows the Coulomb energy coefficient defined as 
ac = EcA^^^Z~'^ with different models as a function of {N — Z)/2 for Pb isotopes. Here, 
Ec denotes the calculated Coulomb energy of a nucleus with a certain mass model. The 
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|N-Z|/2 (N-Z)/2 



FIG. 9: (Color online) (a) Deviations of the calculated b coefficients with DZIO, ETF2 and WS3 
models from bf,t = 0.710^2/3 - 0.946 [16], as a function of isospin. (b) Coulomb energy coefficient 
Uc = EqA^^^Z"'^ obtained with different models as a function of (A^ — Z)/2. 

open circles and the solid line denote the values of Oc with DZIO and WS3, respectively. 
The squares and the crosses denote the results from the HFBTHO [37] and the ETF2 [28] 
approach with the Skyrme force SkP [29], respectively. Similar to Fig. 9(a), the trends 
are quite different with different models. The coefficient Oc in this work is a constant. In 
DZ mass models, the isospin effect in the Coulomb energy calculation is involved by the 
isospin-dependent charge radius Rc oc A^/'^{1 — 0.25/^), which causes the reduction of the 
effective coefficient Oc for neutron- rich nuclei (see the open circles). On the contrary, the 
values of Oc obviously increase with increasing neutron number in the Skyrme energy- density 
functional calculations. It is because for extremely proton-rich nuclei the surface diffuseness 
of protons increases sharply to reduce the huge Coulomb repulsion and causes the reduction 
of Qc consequently. Comparing Fig. 9(a) to (b), one finds that the change of h^i — h with the 
isospin T is roughly consistent with that of the Coulomb energy coefficient ac- The Coulomb 
energy coefficient might be used as a sensitive quantity to study the mass models and the 
isobaric multiplet mass equation. 

In addition, for exploring the correlation among the IMME, the GK relation and the ac- 
curacy of the nuclear mass models, we systematically study the rms deviations of calculated 
h coefficients with Eq.(16) and the nuclear masses by 17 nuclear mass models, including 
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FIG. 10: Values of a{GK) /a{b) as a function of a{GK)/a(M) from 17 different models (solid 
circles). Dashed lines are to guide the eyes. 

different versions of liquid drop models [8, 9, 38, 39], the Dufio-Zuker mass models (DZIO 
and DZ31 [14]), the ETF2 approach with different Skyrme forces, the WS, WS* and WS3 
models. We use a{b) for describing the rms deviation of the calculated b coefficients by a 
certain mass model from bat for the 179 neutron drip line nuclei. The rms deviations a{b) 
obtained from ETF2, DZIO, DZ31 [14] and WS3 are 2734, 870, 630 and 449 keV, respec- 
tively. The values of a{M) and cr(GK) from some mass models are presented in Table II 
previously. In Fig. 10, we show the values of a{GK)/a{b) as a function of a{GK)/a{M). 
Each of solid circles represents the calculated result from a nuclear mass model and all of 
circles are located in between two dashed lines. A general tendency of a{GK)/a{b) increas- 
ing with a{GK)/a{M) is clearly shown, which exhibits a strong correlation between the 
quadratic form of the IMME and the accuracy of the mass model when the GK relations are 
reproduced reasonably well. In other words, reducing the deviation from the IMME could 
be a way to improve the accuracy of mass predictions if the cr(GK) ~ 100 keV. Fulfilling 
both the IMME and the GK relations seems to be two necessary conditions to improve the 
reliability of mass calculations for known nuclei and the extremely neutron-rich nuclei. 
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IV. SUMMARY 



We proposed an extension of our earlier global nuclear mass model which significantly 
improves its ability to describe nuclear masses across the periodic table and puts it at least 
on a par with the very best empirical mass models on the market. With some additional 
correction terms, including mirror nuclei effect due to the isospin symmetry, Wigner-like ef- 
fect from the symmetry between valance-neutron and valance-proton, corrections for pairing 
effects and those from triaxial deformation of nuclei, the rms deviations from 2149 measured 
masses and 1988 neutron separation energies are significantly reduced to 336 and 286 keV, 
respectively. As a test of the extrapolation of the mass model, the a-decay energies of 46 
super-heavy nuclei have been systematically studied. The rms deviation with the proposed 
model reaches 248 keV, much smaller than the result 936 keV from the DZ28 model. Fur- 
thermore, through studying the rms deviations to the masses of nuclei approaching the drip 
lines and the predictive power of the mass formulas fitted to the AME1995 data, we find that 
the results from the proposed model are satisfactory. In addition, with a systematic study 
of 17 global nuclear mass models, we find that the quadratic form of the IMME is closely 
related to the accuracy of nuclear mass calculations when the Garvey-Kelson relations are re- 
produced reasonably well. Fulfilling both the IMME and the Garvey-Kelson relations seems 
to be two necessary conditions to improve the quality of the model prediction. Furthermore, 
the a-decay energies of super-heavy nuclei should be used as an additional constraint on the 
global nuclear mass models. 
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